Decoherence, pointer engineering and quantum state protection 
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We present a proposal for protecting states against decoherence, based on the engineering of pointer 
states. We apply this procedure to the vibrational motion of a trapped ion, and show how to protect 
qubits, squeezed states, approximate phase eigenstates and superpositions of coherent states. 
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It is well known that the interaction of a quantum sys- 
tem with its surrounding environment may lead to quan- 
tum entanglement between system and environment, and 
to an irreversible loss of information on the system. 
Which set of states is less sensitive to entanglement de- 
pends on the concrete form of the interaction Hamilto- 
nian between system and environment |l| . In the course 
of the interaction, the reduced density operator of the 
system becomes rapidly diagonal in this basis, transform- 
ing any initial superposition of these states into a statisti- 
cal mixture. On the other hand, if the system is initially 
in a pointer state, it will remain in a pure state during 
its time development. The decoherence process by which 
coherent superpositions of pointer states get transformed 
into statistical mixtures is at the heart of the quantum 
theory of measurement [^ , and plays an essential role in 
the classical limit of quantum mechanics p| . 

Fighting decoherence has become a major chalenge in 
the last years, motivated by recent progress in the the- 
ory of quantum information processing, which relies on 
the possibility of preserving quantum coherence [Q|5| . It 
is also of interest to high-precision frequency measure- 
ments in ion traps H. Several strategies have been de- 
vised. They include quantum error correction schemes 
H , feedback implementations MM , the realization of q- 
bits in symmetric subspaces decoupled from the environ- 
ment Q , and dynamical decoupling techniques [Ij] . 

In linear ion traps, by far the most important deco- 
herence effect is the one associated with the motional 
state [ p^p^ . In the present paper, we show that de- 
coherence in the vibrational motion of a trapped ion 
can be suppressed by generating, through the techniques 
of "reservoir engineering" pd| , artificial reservoirs asso- 
ciated with properly chosen pointer observables, which 
have the states to be preserved as their eigenstates, and 
which dominate over other dissipation processes. We ex- 
emplify this procedure by showing how to protect from 
decoherence several kinds of non-classical states. 

Under the hypotheses of Markovian dynamics, com- 
plete positivity and initial decoupling between system 
and bath |la,n6[, a master equation describing the re- 
duced dynamics of a system interacting with its enviro- 
ment can be written in the Lindblad form 



where p is the reduced density operator of the system in 
the interaction picture, and we have neglected the unitary 
evolution term —{i/?i)[H, p\. The operators Ci are closely 
related to the system operators present in the interaction 
Hamiltonian and 7^ measures the strength of the system- 
enviroment coupling. In this case, the pointer basis is 
given by the set of all the eigenstates of the operator 
Ci. If all Ci are Hermitian, their eigenstates are steady 
states of the master equation (|l|). If instead the c^'s are 
not Hermitian, the states of the pointer basis will remain 
pure, but will not necessarily be steady states of (|^). 
Note that steady states of P) are not affected by the 
environment ||lO| . Our strategy for the protection of a 
specific quantum state against the enviroment consists in 
engineering, via adequate external driving of the system, 
a system-enviroment coupling, so that the net effect is 
to add to the master equation (nh an extra term, thus 
getting 
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The operator d is chosen so that the state one wants 
to protect is the only steady state of Eq. (||) without 
the environment term Cp. For Pong ^ 7, the steady 
state of the new master equation (g) will be very close to 
the state to be protected (if the state is not unique, the 
term Cp could still induce transitions between the steady 
states). Besides, any state of the system will decay into 
the state chosen to be protected. Therefore, this is also a 
procedure for preparing quantum states in the presence 
of decoherence. 

In the following we apply this method to the one- 
dimensional motion of the center-of-mass of an ion con- 
fined in an electromagnetic trap. The reservoir engineer- 
ing process will be implemented by letting the ion in- 
teract with several laser beams of adequate frequencies 
and intensities, which are quasi- resonant to an electronic 
transition |1) <-» |2) of frequency (^21. 

For our purposes it is important to consider the ion 
to be in the regime of resolved sidebands, given by i^ S> 
P,r2„, where P is the electronic energy decay rate, v is 
the ion vibrational frequency, and ri„ is the (complex) 
Rabi frequency corresponding to laser n, tuned to the 
k-th red vibrational sideband of the ion. Under these 
conditions the interaction term corresponding to laser n, 
can be described, in the interaction picture, by |1^ 
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where g is taken to be real and 
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Here, the operators a and A21 are the annihilation op- 
erator of a quantum of the ionic vibrational motion and 
the electronic flip operator, respectively. The quantity 
rin = •\/?i(k„ • uY /2Mv is the Lamb-Dicke parameter 
with respect to the direction of vibration, fixed by the 
unit vector u. M is the ion's mass and k„ is the wave 
vector of laser n. It is assumed that 77 ^ 1 for any direc- 
tion orthogonal to u, as it is the case in linear traps. 

In the Born-Markov limit, the time evolution of the 
vibronic density operator p for the direction u (say x) is 
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where the second term corresponds to spontaneous emis- 
sion with energy relaxation rate F, and 
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accounts for changes of the vibrational energy along the 
X direction due to spontaneous emission with angular dis- 
tribution W{s). The last term of (g) describes the cou- 
pling of the center-of-mass motion to the environment, 
and has the general form (m). However, the precise form 
of this dissipation term is not important for our purposes. 
The matrix elements of (||) with respect to the elec- 
tronic basis yield the equations 
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We assume now that the decay rate F is by far the 
largest rate in the system. Under this condition, one can 
eliminate pi2 adiabatically. Since £ ex 7, one gets: 

P12 = -{2ig/T){d)p22 - PiiJt) [1 -f O (7/F)] . (10) 

Replacing (^Ol) into (^) and (pi) , and adding up these two 
equations, we get, since the reduced density operator for 
the vibrational mode is given by pi, = pn + P22, and 
neglecting the correction proportional to 7/r in (M): 
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Under the conditions assumed here, the matrix ele- 
ments of P22 are much smaller than those of pn. Indeed, 
replacing (^ in (||), and eliminating P22 from (^) adia- 
batically, one gets P22 ~ O [(g/F)^] pn. We can there- 
fore safely neglect the terms dependent of P22 inside the 
brackets in (^_l|) and at the same time replace pn by p^. 
We verified numerically that these are indeed excellent 
approximations. We get then, finally: 
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We will base our considerations on this equation. The 
first term on the r.h.s. has the form (m. This is the "en- 
gineered reservoir," with a decay constant Feng = 4:g'^/^- 

Neglecting terms of 0(77"*) in the expansion of the sec- 
ond term on the r.h.s. of eq. ([l2|), one can show that its 
contribution is ~ {ri'^/b){4:g^/r)py, that is, (277^/5) mul- 
tiplied by the engineered-reservoir term. For rj « 0.25, 
this yields a factor ~ 1/40, a small correction, which is 
however fully taken into account in our numerical simu- 
lations. Therefore, the action of the engineered reservoir 
will be the dominant one as long as Feng 3> 7. 

In recent experiments with trapped ions, random fields 
seem to play an important role in the decoherence pro- 
cess [Q. Their effect may also be described by (p. We 
write the random field as E{t) ^ £(+)(i) exp(— wt) + 
£^~'{t) exp{iiy't), and the interaction Hamiltonian in the 
RWA as ||: H = -p[S'^+'){t)a'i + £^-\t)a]. Iter- 
ating the equation of motion for the density operator, 
and using that for steady fields |l8) {£^-'^{t)£'--\t')) = 
(£(+'(<)£(+)(*')) = 0, we get, in the Markovian limit 
(£W(t)£(-)(t')> ^2DS{t^t'): 
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which corresponds to an infinite temperature thermal 
reservoir (by letting the thermal photon number Nt — > 
cxD, and at the same time the dissipation rate 7 — > 0, 
so that jNt remains constant). Both random fields and 
thermal reservoirs will be considered in our simulations. 

In order to protect a state {ip), we look for d such that 
d\ip) — X\ip) with A = 0, and make sure that \ip) is in fact 
the only steady state of (0) without Cp. Since Feng ^ 7, 
this yields a good approximation of the corresponding 
steady state of (0). 

As a first example, we consider the protection (and 
generation) of the class of states \ip) = X]n=o '^n\n), where 
1 71) is an energy eigenstate of the vibrational motion of 
the trapped ion and c„ ^ 0. It is easy to see that the 
operator d — g(n) a + h{n) has \ip) as its only eigenstate 
with eigenvalue A = 0, provided the eingenvalues of giji) 
and 11(77,) fulfill the constraints g{m) — — / ™ , Cm/cm+i 
{m — 0,- ■ ■ , N— 1) and N is the first zero of h(77T,). For 
this case, Jt/)) is the only steady state of (0). 



Inspection of Eqs. (||) and (Q) shows that the operator 
g(n) a can be engineered by driving the ion with N laser 
fields, tuned to the first vibrational sideband of the ion. 
The values of the Rabi frequencies fin of the N lasers 
are given by the following system of N linear equations 
(m = 0, •••,7V-1): 
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where the Lamb-Dicke parameters ??„ depend on the ori- 
entation of the laser fields with respect to the x-direction. 
The operator h(n) is constructed by driving the ion with 
two laser fields resonant with the electronic transition, 
one of them propagating orthogonal to the x-axis (say 
y-direction, with rjy <^ 1). The Rabi frequencies of these 
two lasers are related by fly = —Q,xLi^{tj^), where Ln{x) 
is a Laguerre polynomial of order iV (rjx should not be 
too large for the first zero of /i(to) to occur at m = N). 
An important representative of the class of states pre- 
sented above is the "qubit" state \tp) = co|0)-|-ci|l). The 
discussion above implies that this state can be generated 
and protected against the action of an external reservoir 
with just three lasers, with Rabi frequencies satisfying 
the following condition: 
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In this case, Fong = TfVL\/T in Eq. (||). In order for 
the corresponding reservoir to win over the environment 
reservoir, one needs rf'Vli /T 3> 7, but at the same time 
r, J7i <C z/ and T 3> "qfli- This requirements are satisfied 
if r « 4 MHz, f7i w 2 MHz, 77 = 0.2, ;/ w 20 - 30 MHz, 
as long as 7 < rfVli'^/T w 40 kHz. Fig. |i|a displays 
the fidelity ¥{t) = Tr{/5(0)p(t)}, with the ion initially in 
the vibrational state jV') — jdO) + |1)) [e^U our numerical 
simulations are obtained from eq. (pf), and we always as- 
sume the ion to be initially in the electronic ground state] . 
Both a thermal and a random field reservoir have been 
considered. As can be seen, the system rapidly reaches a 
steady state with fidelity very close to unity (77 = 0.2). 

One should remark that feedback procedures Q do 
not protect states involving superpositions of |0) and |1), 
since the loss of one photon by the state |1) completely 
erases any phase information about the original state. 
Our procedure works however very well in this case. 

Other interesting example is the approximate phase 
state H IV) = (l/VFTT) ELo e™-*!"), which can be 
generated and protected by A^-|-2 lasers. Fig. |^b displays 
the time evolution of the fidelity F(t), for the approx- 
imate phase state |i/') = (1/2) ELo ^'""^1") (^'^ i^ ^^le 
range 0.1 - 0.2). 

A class of states which is specially fragile against the 
action of decoherence is the one formed by mesoscopic su- 
perpositions of coherent states. Under action of an exter- 
nal reservoir, these states decay to a mixture of coherent 



states in a extremely short time, inversely proportional 
to the distance between the two states in phase space [^ . 
Our technique can also be applied to the Schrodinger- 
cat-like state (20| |<?!)+) = (|a) +i\- a))/\/2. Since \(f)+) 
has no "hole" in its number distribution (which coincides 
with that for the coherent state \a)), it can be approxi- 
mated by one of the states 1-0) discussed above. Conse- 
quently, one could, with the use of TV -)- 2 lasers, generate 
and protect the state li/)) — X]n=o "^"l"")' ^'^^^ the first 
N coefficients c„ equal to the corresponding coefficients 
of 10+) . One should notice, however, that it is possible 
in this case to find directly a Lindblad operator d, which 
has the state |0+) as its only eigenstate with zero eigen- 
value: d = TdT'^ = e""a -I- ia. Here T is the unitary 
operator exp[i7rn(ri — l)/2] exp(aa^ — a*a), which yields 
I (/)-(-) when applied to the vacuum. For this choice of d, 
we plot in Fig. Ic the fidelity F{t), for the initial state 
\4'+) iv — 0.2). An open problem is how to engineer this 
operator with a finite number of laser beams. 

Finally, we describe the protection of a squeezed state. 
We set d = a + xd\ where x = tanhr and r is the 
squeezing factor. The corresponding setup consists of two 
lasers along the direction of squeezing, resonant with the 
first red (laser 1) and the first blue (laser 2) sidebands, 
and with Rabi frequencies satisfying fl2/fli = x p4[ - 
The numerical simulation is shown for r — 0.6 in Fig. Id, 
for a realistic set of parameters (77 = 0.05). Higher values 
of squeezing render our method less effective, since the 
presence of higher photon numbers would lead to stronger 
dissipation by the "natural" reservoir. 

In conclusion, we have suggested a method for protect- 
ing quantum states of the vibrational motion of a trapped 
ion against decoherence by generating artificial reservoirs 
which have the states to be protected as pointer states. 
More general pointer states can be generated by applying 
unitary transformations to the states and operators dis- 
cussed here. Indeed, the transformed states would still be 
the sole steady solutions of the master equation with the 
transformed operators. This is precisely the mechanism 
which leads to the protection of the states |0+) above, 
and also of squeezed states, since they are related by 
unitary transformations to the vacuum, which is the only 
steady state for a zero-temperature reservoir (for which 
d — a). As for possible sources of error in experimen- 
tal implementations of our proposal, one should remark 
that our method is sensitive only to relative intensity and 
phase fluctuations, since state selection is determined by 
ratios of Rabi frequencies. 
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FIG. 1. Time evolution of the fidelity F{t) for several initial 
nonclassical states, in the presence of a thermal a^ = 3; and 
(d) vacuum squezed state with r — 0.6. 
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